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ABSTRACT 

We present MMT/Megacam imaging of the Leo IV dwarf galaxy in order to investigate its structure 
and star formation history, and to search for signs of association with the recently discovered Leo V 
satellite. Based on parameterized fits, we find that Leo IV is round, with e < 0.23 (at the 68% 
confidence limit) and a half-light radius of ~ 130 pc. Additionally, we perform a thorough search 
for extended structures in the plane of the sky and along the line of sight. We derive our surface 
brightness detection limit by implanting fake structures into our catalog with stellar populations 
identical to that of Leo IV. We show that we are sensitive to stream-like structures with surface 
brightness fi r < 29.6 mag arcsec -2 , and at this limit, we find no stellar bridge between Leo IV (out to 
a radius of ^0.5 kpc) and the recently discovered, nearby satellite Leo V. Using the color magnitude 
fitting package StarFISH, we determine that Leo IV is consistent with a single age (~14 Gyr), single 
metallicity ([Fe/H] ~ —2.3) stellar population, although we can not rule out a significant spread in 
these value. We derive a luminosity of My — —5.5 ± 0.3. Studying both the spatial distribution 
and frequency of Leo IV's 'blue plume' stars reveals evidence for a young (~2 Gyr) stellar population 
which makes up ^2% of its stellar mass. This sprinkling of star formation, only detectable in this deep 
study, highlights the need for further imaging of the new Milky Way satellites along with theoretical 
work on the expected, detailed properties of these possible 'reionization fossils'. 
Subject headings: 



1. INTRODUCTION 

Since 2005, 14 satellite co mpanions to the Milky Way 
have been discovered (see IWillmanl (|2010D and refer- 
ences therein). Despite the fact that many of these ob- 
jects are less luminous than a typical globular cluster 
(—1.5 < My < —8.6), these 14 objects have a range 
of properties that encompass the most extreme of any 
galaxies, i ncluding: the highest inferred dark matte r con- 
tent (e.g. iSimon fc QenaH2007t iGeha et all [2009h . the 
lowest [Fe/H] content ( Kirbv et al.ll2008l), unusually ellip- 
tical morphologi es (e.g. Hercules: IColeman "eTaLl [20071: 
ISand et alJl2009h . and in some cases evidence for severe 
tidal disturbance ( e.g. Ursa Major II; IZucker et al.ll200d 
IMunoz et al.l[20 09). The varied properties of these lowest 
luminosity galaxies are valuable probes for understand- 
ing the physics of dark matter and galaxy formation on 
the smallest scales. 

Of the newly discovered Milky Way (MW) satellites, 
Leo IV {My = —5.5, ~ 130 pc) is among the least 
studied, despite several signs that it is an intriguing 
object. Leo IV appears to be dominated by an ol d 
and metal-pool stellar population (|de Jong et al.l [2~008b). 
However it also has an apparently complex color magni- 
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tude diagram (CMD), with a 'thick' red giant branch, 
possibly caused by either multiple st ellar populations 
and / or depth along the line of sight (jBelokurov et al.1 
2007). Leo IV also may have a very extended stellar dis- 
tribution, despite its apparently round (e = 0.22l ; 22 ) 

and compact (r/j = 2.5^7 arcmin; iMartin et aLll2008bl ) 
morpholog y. A search for variab le stars was recently per- 
formed by iMoretti et al.l ([2009D . who used the average 
magnitude of three RR Lyrae stars to find a distance 
modulus of (m—M) — 20.94±0.07 mag, corresponding to 
154±5 kpc. Interestingly, one of the three RR Lyrae vari- 
ables lies at a projected radius of ~10 arcmin, roughly 
three times the half light radius, leading to the suggestion 
that Leo IV may actually possess a 'deformed morphol- 
ogy'. 

Based on Keck/DEIMOS spectroscopy of 18 member 
stars, Leo IV has one of the smallest velocity dispersions 
of any of the new MW sa tellites, with a = 3.3 ± 1.7 
km/s ([Simon fe Gehall2007l). A m etallicity study of 12 
of these spectra ([Kirbv et al.ll2008l ) showed Leo IV to be 
extremely metal poor, with ([Fe/H]) = —2.58 with an 
intrinsic scatter of <J\Fe/H] = 0.75 - the highest disper- 
sion among the new dwarfs. 

Recently, the MW satellite Leo V (My 4.3 ± 0.3) 

has been discovered, separated by only ~2.8 degrees on 
the s ky and ~40 km/s from Leo IV (Belo kurov et al.1 
|2008| ) . With a Leo V dist ance of ~180 kpc, this cl ose sep- 
aration in phase space led lBelokurov et all (|2008j) to sug- 
gest that the Leo IV/Leo V system may be physically as- 
sociat ed. This argument was bolstered bv lWalker et al.l 
(2009), who spectroscopically identified two possible 
Leo V members 13 arcminutes from the satellite's center 
(Leo V's rh is ~0.8 arcminutes) along the line connecting 
Leo IV and Leo V, suggesting that Leo V is losing mass. 
A recent analysis bv lde Jong et al.l ([20101 ) of two 1 square 
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degree fields situated between Leo IV and Leo V shows 
tentative evidence for a stellar 'bridge' between the two 
systems with a surface brightness of ~32 mag arcsec -2 . 

Motivated by all the above, we obtained deep photom- 
etry of Leo IV with Megacam on the MMT. In this paper, 
we use these data to present a detailed analysis of both 
the structure and SFH of Leo IV. We also search for any 
signs of disturbance in Leo IV which may hint at a past 
interaction with the recently discovered, nearby Leo V. 
In § [5] we describe the observations, data reduction and 
photometry. We also present our final catalog of Leo IV 
stars. In §|3]we derive the basic structural properties of 
Leo IV, and search for signs of extended structure. We 
quantitatively assess the stellar population of Leo IV in 
§ H] using both CMD-fitting software and an analysis of 
its blue plume population. We discuss and conclude in 

§13 

2. OBSERVATIONS AND DATA REDUCTION 

We observed Leo IV on April 21 2009 (UT) with Mega- 
cam (|McLeod et alj2000h on the MMT. MMT/Me gacam 
has 36 CCDs, each with 2048 x 4608 pixels at O'.'08/pixel 
(which were binned 2 x 2), for a total field-of-view (FOV) 
of ^24'x24'. We obtained 5 250s dithered exposures in 
g and in r in clear conditions with Of.' 9 seeing. We re- 
duced the data using the Megacam pipeline developed 
at the Harvard-Smithsonian Center for Astrophysics by 
M. Conroy, J. Roll and B. McLeod, which is based in 
part on M. Ashby's Megacam Reduction Guide 8 . The 
pipeline includes standard image reduction tasks such 
as bias subtraction, flatficlding and cosmic ray removal. 
Precise astrometric solutions for each science exposure 
were derived using t he Sloan Digital Sky Survey Dat a Re- 
lease 6 (SDSS-DR6 lAdelman-McCarthv et alJl2006l ) and 
the constituent images were then resampled onto a com- 
mon grid (using the lanczos3 interpolation function) and 
combined with SWarp 9 using the weighted average. 

Stellar photometry was determined on the final im- 
age stack using a nearly identical methodology as 
ISand et "all (|2009l ) with the comma nd line version of the 
DAOPHOTII/Allstar package ([Stetson| Gill) ■ We 
allowed for a quadratically varying PSF across the field 
when determining our model PSF and ran Allstar in 
two passes - once on the final stacked image and then 
again on the image with the first round's stars sub- 
tracted, in order to recover fainter sources. We culled 
our Allstar catalogs of outliers in \ 2 versus magni- 
tude, magnitude error versus magnitude and sharpness 
versus magnitude space to remove objects that were not 
point sources. In general, these cuts varied as a func- 
tion of magnitude. We positionally matched our g and r 
band source catalogs with a maximum match radius of 
0'.'5, only keeping those point sources detected in both 
bands in our final catalog. 

Instrumental magnitudes are put onto a standard pho- 
tometric system using stars in common with SDSS-DR6. 
We used all SDSS stars within the field of view with 
17.5 < r < 21.5 and -0.25 < g - r < 1.5 to perform the 
photometric calibration and simultaneously fit for a ze- 
ropoint and a linear color term in g — r. The linear color 
term slope was 0.11 in (r — tmmt) versus (g — r) and 

8 http:/ /www. cfa.harvard.edu/~mashby/megacam/mcgacam_frE 

9 http:/ /astromatic. iap.fr/software/swarp/ 



0.086 in (g — jmmt) versus (g — r), consist ent with that 
found in the MMT study of Bootes II by IWalsh et all 
(2008). Slight residual zeropoint gradients across the 
field of view were fit to a quadratic function and cor- 
rected for (see also Saha et al. 2010), resulting in a final 
overall scatter about the best-fit zeropoint of 5g ~0.05 
and Sr ~0.05 mag. 

To calculate our photometric errors and completeness 
as a function of magnitude and color, we performed a se- 
ries of artificia l star tests using a technique nearly iden- 
tical to that of ISand et al.l (|2009h . Briefly, artificial stars 
were placed into our Leo IV images on a regular grid 
with spacing between ten and twenty times the image 
full width at half maximum with the DAOPHOT rou- 
tine ADDSTAR. In all, ten iterations were performed 
on our Leo IV field for a total of ~ 300000 implanted ar- 
tificial stars. The r magnitude for a given artificial star 
was drawn randomly from 18 to 29 mag, with an ex- 
ponentially increasing probability toward fainter magni- 
tudes. The g — r color is then randomly assigned over the 
range -0.5 to 1.5 with equal probability. These artificial 
star frames were then run through the same photometry 
pipeline as the unaltered science frames, applying the 
same x 2 > sharpness and magnitude-error cuts. For ref- 
erence, we are 50% (90%) complete in g at ^25.3 (23.6) 
mag and in r at 24.8 (23.3) mag. When necessary, such 
as for calculating the SFH of Leo IV in § 0] the complete- 
ness as a function of both magnitude and color is taken 
into account. 

2.1. The Color Magnitude Diagram and Final Leo IV 

Catalog 

We present the CMD of Leo IV in Figure [TJ Plotted in 
the left panel are all stars within the half-light radius (as 
determined in § 13. ip , while in the right panel we present 
a Hess diagram of the same field with a scaled back- 
ground subtracted using stars located outside a radius of 
12 arcminutes. 

In both panels we highlight the possible stellar popu- 
lations of Leo IV. In th e right panel we plo t three the- 
oretical isochrones from iGirardi et ail (|2004| ). The solid 
and dashed lines are old 14 Gyr isochrones with [Fe/H] 
of -2.3 and -1.7, respectively. We adjust these isochrones 
to have m — M=20.94, as found in the RR Lyrae study 
of Leo IV bv iMoretti et al.l (|2009fi . which fit the ridge- 
line well. With the dotted line, we also plot a 1.6 Gyr 
isochrone with [Fe/H]=-1.3. This isochrone agrees rela- 
tively well with the 'blue plume' stars which are evident 
in Leo IV's CMD. 

We mark several regions in the left panel of Figure Q] 
for quantitative study in later sections. The solid box 
denotes the blue horizontal branch (BHB) stars, which 
are clearly defined in our CMD. The dashed and dotted 
regions are two different selection regions for blue plume 
stars. The dashed is the total blue plume population, 
although much of this region may be plagued by fore- 
ground stars or unresolved galaxies. We thus will also 
utilize the stars within the dotted box as a relatively 
contamination-free tracer of the blue plume star popula- 
tion in later sections. An open question is whether or not 
the blue plume stars are young stars, as is plausible based 
on the CMD, or are blue stragglers. We quantitatively 
Hllgess the blue plume population in § 14.21 

Given that the CMD of Leo IV is visually in ex- 
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cellent agreement w ith the distance measurement of 
iMoretti et al.l ()2009f ) . and the possible presence of mul- 
tiple stellar populations, we do not attempt a CMD- 
fittin g method for m easuring the distance to Leo IV 
(e.g- ISand et al.ll2009f ). and will adopt (to - M) = 20.94 
throughout this work. We will vary this quantity when 
necessary to determine how sensitive our results are to 
this assumption. 

We present our full Leo IV catalog in Table Q] which 
includes our g and r band magnitudes (uncorrected for 
extinction) with their uncertainty, along with the Galac- 
tic ex tinction values derived for each star (Schlegel et aL 
1998). We also note whether the star was taken from the 
SDSS catalog rather than our MMT data, as was done for 
objects near the saturation limit of our Megacam data. 
Unless stated otherwise, all magnitudes reported in the 
remainder of this paper will be extinction corrected. 

3. LEO IV STRUCTURAL PROPERTIES 

We split our analysis of the structural properties of 
Leo IV into two components. First, we fit parameterized 
models to the surface density profile of Leo IV. Following 
this, we search for signs of extended structure in Leo IV, 
especially in light of its proximity to Leo V. 

3.1. Parameterized Fit 

It is common to fit the surface de nsity profile of both 
globu lar clusters an d dSphs to King (|Kindll966l ). Plum- 
mcr (PI ummeriri911| ). and exponential profiles. This is 
an important task both to facilitate comparisons with 
other observational studies and f or understanding th e 
MW satellites as a population (e.g.. lMartin et al. 2008b). 
To do this, we use the CMD selection region shown in 
Figure [5] for isolating likely Leo IV members. This CMD 
selection box was determined by first taking a M92 ridge- 
line at (to — M)=20 .94, the distance to Leo IV found by 
IMoretti et al.l (|2009l ). and placing two bordering selec- 
tion boundaries a minimum of 0.1 mag on either side in 
the g — r color direction. These selection regions are in- 
creased to match the typical g — r color uncertainty at a 
given r magnitude (as determined with our artificial star 
tests) when that number exceeds 0.1 mag. A magnitude 
limit of r=24.8 mag was applied to correspond to our 
50% completeness limit. At the present we focus on Leo 
IV ridge line stars, but we discuss the spatial properties 
of the blue plume and horizontal branch stars in § 14.21 

We fit the three parameterized density profiles to the 
stellar distribution of Leo IV: 
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where rp and a are the scale lengths for the Plummer 
and exponential profiles and r c and r t are the King core 
and tidal radii, respectively. For the Plummer profile, rp 
equals the half-light radius r^, while for the exponential 
profile 77, ~ 1.668a. We simultaneously fit a background 



surface density, £t>, while fitting the Plummer and expo- 
nential profiles. For the King profile, there is a degener- 
acy between the tidal radius and the background surface 
density. We thus fix the background value to the average 
of that found for the Plum mer and exponentia l profiles 
for our King profile fits (e.g. IWalsh et al.ll2008t ). 

We use a maximum likelihood (ML) technique for 
co nstraining str u ctural parameters similar to that 
of IMartin et al.l (|2008bf ). which we have refined in 
ISand et al.l ( 20091) . and further refined in the current 
work. We point the reader to those works for further 
details concerning the expression of the likelihood func- 
tion. Including the central position, ao and Sq, posi- 
tion angle (6), and ellipticity (e) both the exponential 
and Plummer profiles have the same free parameters - 
(ao,5Q,0,e,rhaif,^b), while the King profile free param- 
eters are (oto,5o,9,e,r c ,rt). Uncertainties on structural 
parameters are determined through 1000 bootstrap re- 
samples, from which a standard deviation is calculated. 
We have tested the robustness of our algorithm for dwarf 
galaxies with roughly the same number of stars as Leo 
IV in an Appendix, and will use these tests to inform our 
results in what follows. 

Our results are presented in Table [5] We show our 
best fit stellar profiles in Figure [3l Although the plot- 
ted stellar profiles are not fit to the plotted binned data 
points, they do show excellent agreement. We note that 
the apparent slight overdensity at R ^8 arcmin above 
the derived parameterized fits does not correspond to 
any single feature, as can be seen from the smoothed 
map of Leo IV that we present in § l3.2l and Figure HI but 
is likely just the result of several fluctuations at roughly 
the same radius. Interestingly, Leo IV appears to be par- 
ticularly round, at least according to the parameterized 
model fit to the data. In fact, our ML-derived ellipticity 
is consistent with zero (see also the Appendix), allowing 
us to only place an upper limit of e < 0.23 (at the 68% 
confidence limit). Given this low ellipticity, we can not 
place any meaningful constraints on the position angle of 
Leo IV, as both the tests in the Appendix indicate and 
the bootstrap resamples of our Leo IV data reaffirm. On 
a separate note, the tidal radius for the King profile fit, 
with a value of r t — 18.55', is larger than the Megacam 
FOV. Thus, this value should be taken with caution. 

It is useful to compare our parameterized fit results 
with simila r work in the literatu re. Using the original 
SDSS data, IMartin et all (|2008bf ) fit Leo IV with an ex- 
ponential profile using a ML technique similar to the one 
utilized in the current work, and found re sults within 1- 
a of t hose presented here. More recently, Ide Jong et al.l 
(2010) used deeper data from the 3.5 m Calar Alto tele- 
scope around both Leo IV and Leo V, again applying a 
ML algorithm to measure their structure. In this case, 
the authors find r% = 4.6+^7 an d a well measured ellip- 
ticity of e — 0.49 ± 0.11, both statistically inconsistent 
with the resu lts presented in the c urrent work. As com- 
mented on by Ide Jong et al.l (|2010[) . this is likely because 
of the different stellar populations probed - while the 
present work uses mostly main sequence and subgiant 
stars t o determine the structure of Leo IV, Ide Jong et al.l 
(2010) use mostly brighter stars and objects on the blue 
horizontal branch. Also, Leo IV lies on the edge of 
one of their pointings, which could have biased their re- 
sults. Resolution of this discrepancy will require addi- 
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tional measurements. 

3.2. Extended Structure Search 

We now search for signs of tidal disturbance and other 
anomalies - which would not be picked up by our pa- 
rameterized fits in § 13.11 - based on the morphology of 
Leo IV's isodensity contours. We do this with an eye 
towards determining if there is a current physical con- 
nection between Leo IV and Leo V. 

Our basic approach is similar to that of ISand et"al"1 
(2009) . We include stars within the same CMD selection 
box used for our parameterized fits (Figure [2]) , placed 
those stars in 10" x 10" bins and spatially smoothed these 
pixels with three different Gaussians of <r=0.5,1.0 and 1.5 
arcminutes. The background level of stars in the CMD 
selection box, and its variance, was determined via the 
MMM routine in IDL 10 . To avoid the bulk of Leo IV, 
these statistics were determined in two boxes of size 20' x 
3' situated 9.5 arcminutes North and South of Leo IV. 
We found that our smoothed maps were unaffected if 
only one of these boxes was used, or if we varied their 
sizes. We present our smoothed maps in Figure HI with 
the contours representing regions that are 3, 4, 5, 6, 7, 
10 and 15 standard deviations above the background. 
We focus on the <r=0.5 arcminute map, since it contains 
the most detail without loss of potential Leo IV features. 
The cr=1.0 arcminute smoothing scale will be useful in 
§ 13.2.11 as we explore our sensitivity to stellar streams. 

Outside of the main body of Leo IV, there appears to 
be only a handful of compact, 3, 4 and 5 a overdensi- 
ties. How significant are these overdensities, given the 
binning, smoothing and number of pixels that went into 
the making of Figure HP To gauge their significance, 
we take our input photometric catalog, randomize the 
star positions and remake our smoothed maps (with the 
0.5 arcmin Gaussian) for several realizations (Figure [5]). 
While some of these maps have several 3 and 4 a over- 
densities, while others have fewer, we find that the distri- 
bution of pixel values maintains a Poisson distribution. 
We thus conclude that the majority of features outside 
the main body of Leo IV in Figure [4] are likely just noise 
- with the possible exception of the 5-cr overdensities at 
positions ( 6',- 12') and ( 8', 10') with re- 
spect to Leo IV. Background-subtracted Hess diagrams 
of these two regions were made from our Leo IV catalog, 
but they do not yield CMDs that are consistent with 
Leo IV's stellar population (see § 13 . 2. If) . Thus, our ob- 
servations yield no strong evidence for substructure in 
the vicinity of Leo IV. 

The main body of Leo IV itself has some interesting 
features. There is a hint of an elongation or disturbance 
in the core of Leo IV, along with two 'tendrils' - one 
directed to the West and the other to the Southwest. 
Again, due to the small number of stars in Leo IV, these 
irregularities may be an effect of small number statistics. 
To evaluate the significan ce of the morpholo gy in Fig- 
ure HI we follow the path of I Walsh et"al"1 |2008| ) and their 
evaluation of morphological irregularities in Bootes II. 
We bootstrap resample the Leo IV stars and replot our 
smoothed maps. The results of nine such resamples can 
be seen in Figure HJ While we cannot rule out that the 
tendrils seen in our Leo IV map are genuine, they are 

10 available at http://idlastro.gsfc.nasa.gov/ 



not ubiquitous features in our resampled maps, and so 
we cannot with confidence claim they are real features. 

We finally point out that there appears to be no sign 
of interaction or disturbance in the direction of Leo V, 
which we indicate in the middle panel of Figure |4] A 
Hess diagram of all stars farther than 5 arcminutes from 
the center of Leo IV, and within 1.5 arcminutes of the 
line that connects Leo IV and Leo V, yields a differ- 
enced CMD that is consistent with noise (after proper 
background subtraction). Any such disturbance would 
have to be below our detection threshold, which we de- 
termine in the next section. Finally we point out that 
these smoothed maps are sensitive to structures at the 
distance of both Leo IV and Leo V (~180 kpc; Belokurov 
et al. 2009) given that we are predominantly probing the 
CMD at magnitudes brighter than the subgiant branch. 

3.2.1. Inserting Artificial Remnants 

In order to assess our surface brightness limit, and 
our sensitivity to structures of different sizes and mor- 
phologies, we insert fake 'nuggets' and 'streams' into our 
Leo IV catalog with stellar populatio ns drawn from on e 
consistent with that of Leo IV. As in lSand et al.l (|2009h . 
we use the testpop program within the CMD-fitting pack- 
age, StarFISH, to generate our artificial 'Leo IV CMDs 
and then remake our smoothed maps. By varying the 
number of stars (and, by extension, the surface bright- 
ness) in these structures, we can then assess whether 
our adopted search method for extended structure would 
have recovered them, and if so, what the resulting CMD 
would look like. This method of 'observing' and then 
examining these artificial remnants is more informative 
than traditional methods of simply quoting a '3 — a' sur- 
face brightness limit, even if the result is a relatively 
ambiguous detection limit. For d e tails of the procedure, 
we refer the reader to lSand et al.l (|2009l ). and to our pre- 
sentation of Leo IV's SFH in § U 

We inject both 'nuggets' - Leo IV-like stellar popu- 
lations with an exponential profile having rhaif = L0 
arcmin - and 'streams' - with a Gaussian density pro- 
file in the right ascension direction with rr=1.5 arcmin 
and a uniform distribution in the declination direction 
over the Northern half of the field - into our final Leo IV 
photometric catalog. In Figures [7] and we illustrate 
some results from our tests, along with their properties 
in Tabled 

We show an example of a 35 star 'nugget' near what 
we believe is our detection threshold in Figure El While 
this nugget has a peak detection at 6.4-c, it is not par- 
ticularly different in morphology or significance than the 
other random peaks in our Leo IV field. This, however, 
changes when the resulting Hess diagram is examined, 
showing several stars along the red giant branch - some- 
thing that is not apparent in the true overdensities in 
our field. Taking this as our rough detection limit for 
compact remnants of Leo IV, we are sensitive to ob- 
jects with a central surface brightness of /io,r = 27.9 and 
Mo, g = 28.3. 

Turning towards our artificial 'streams', both the 
200 and 300 star streams are easily detectable in our 
smoothed maps (where we have used a 1 arcminute 
smoothing scale to better pick out the thick streams - 
Figure [5]). Further, the resulting Hess CMDs show a 
considerable red giant and BHB presence, and the be- 
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ginnings of the main sequence for the 300 star case. The 
analogous 100 star stream is not convincingly detected. 
We thus suggest that we are reliably sensitive to streams 
with central surface brightness [i r ~ 29.6 and fi g ~ 29.8 
(as measured along the center of the stream) with a ge- 
ometry and morphology roughly similar to that simu- 
lated. 

The recent work of Ide Jong et al.l (|2010D have pre- 
sented tentative evidence for a stream connecting Leo IV 
and Leo V with a surface brightness of ~32 mag arcsec -2 . 
Unfortunately, despite our much deeper data, we are un- 
able to probe down to such faint surface brightness limits 
due to the relatively small area of our single pointing. 

3.3. Structure along the line of sight 

Because previous authors have suggested that the 
'thick' red giant branch in Leo IV may be due to elon- 
gation along the line of sight, it is worth searching for 
signs of such elong ation in the width of the BHB (e.g., 
iKlessen et al.1 12003). To do this, we have created a BHB 
star fiducial using data collected by J. Strader from SDSS 
encompassing 41 horizontal branch stars from M3 and 
M13, corrected for extinction and their relative distances. 
Using these stars, a third-order polynomial was fit with 
the IDL routine robust_POLY_fit, and we used this 
polynomial to define our BHB fiducial. 

For simplicity, we assume that all stars in our Mega- 
cam field within the solid box in Figure [TJ are BHB stars 
belonging to Leo IV. Placing the BHB fiducial at our as- 
sumed distance to Leo IV, (to— M) = 20.94, leads to a vi- 
sually excellent match. However, since we are interested 
in the scatter about this fiducial to put constraints on the 
width of Leo IV, rather than the absolute distance to it, 
we adjust the fiducial BHB sequence so that the average r 
magnitude deviation of the Leo IV BHB stars against the 
fiducial is ~0.0 (this adjustm ent was 0.06 mag, af firming 
the distance measurement of iMoretti et al.1 ([20 09)). The 
resulting root mean square deviation of the Leo IV BHB 
stars about the fiducial is ^0.2 mag, which corresponds 
to a deviation of ~15 kpc at the distance of Leo IV. 
While this limit is comparable to the quoted difference 
in dis tance between Leo IV and Leo V (|Belokurov et al.l 
2008), it is also roughly the measurement limit achievable 
using the spread in BHB mags as a measure of line-of- 
sight depth. First, there are know n RR Lyrae variable s 
among the BHB stars in Leo IV (jMoretti et al.ll2009f) . 
which we have observed at a random phase. Second, 
for a spread in metallicity of <J[Fe/H] = 0-?5 in Leo IV 
(|Kirbv et al.l l20081. one would expec t a natural spread in 
BHB star magnitudes of ~ 0.2 mag (|Sandage fc Cacciaril 
Il990t [Olszews ki et alJ[l996[ ) . We therefore conclude that 
our data, while consistent with no elongation along the 
line of sight, can not put strong constraints on the depth 
of Leo IV. 

4. STELLAR POPULATION 

There is spectroscopic evidence tha t Leo IV has a large 
metallicity spread, and according to iKirbv et al.l (2008) 
it has the largest metallicity spread seen among the new 
dwarfs ({[Fe/ H])=— 2.58 with a spread of cr[ Fe / H ]=0.75). 
Additionally, Belo kurov et al.1 (|2007| ) have hinted that 
Leo IV has a 'thick' red giant branch, indicating a spread 
in metallicity, and this appears to be the case (at least 
superficially) in the CMD presented in the current work 



(Figure [!}. Also intriguing is the population of blue 
plume stars pointed out in Figure Q] and § 12.11 In this 
section, we first perform a CMD-fitting analysis of the 
stellar population of Leo IV with StarFISH, and then go 
on to look at the blue plume population in detail. We 
end the section by combining the structural properties 
found in § 13.11 with the stellar population determined in 
the current section to calculate the total absolute mag- 
nitude of Leo IV. 

4.1. Star Formation History via CMD-fitting 

Here we apply t he C MD-fitting package StarFISH 
([Harris fc Zaritskvl I2001D to our Leo IV photometry 
within the half light radius to determine its SFH and 
metallicity evolution. As discussed in previous works, 
StarFISH uses theoretical isochrones (we use those of 
Girardi et al. 2004, although any may be used) to con- 
struct artificial CMDs with different combinations of dis- 
tance, age, and metallicity. Once convolved with the ob- 
served photometric errors and completeness (using the 
artificial star tests of § [5]) , these theoretical CMDs are 
converted into realistic model CMDs which can be di- 
rectly compared to the data on a pixel-to-pixel basis, af- 
ter binning each into Hess diagrams. We use the Poisson 
statistic of Dolphin] (|2002l ) as our fit statistic. The best 
fitting linear combination of model CMDs is determined 
throug h a modification of the standard AMOEBA algo- 
rithm |PriiIeiaI] [l9H EiEiE&lilitik3[2Q0l|). Sev- 
eral steps are taken to determine the uncertainties in 
StarFISH fits, which are discuss e d in d etail in previ- 
ous wo rk - see lHarris &: Zaritskvl (|2009f ): Ide Jong et al.l 
(|2008aD . 

Our StarFISH analysis is similar to that of I Sand et al.l 
(j2009[ ). We include isochrones with [Fe/H]=-2.3,-1.7 and 
-1.3 and ages between ^10 Myr and ^15 Gyr. Age 
bins of width Alog(t)=0.4 dcx were adopted, except for 
the two oldest age bins centered at ^10 Gyr and ^14 
Gyr, where the binning was Alog(t)=0.3 dex. A 'back- 
ground/foreground' CMD, created by taking all stars 
outside an elliptical radius of 12 arcminutes (with ellip- 
ticity of e=0.05, as in our best-fitting exponential profile, 
see Table [5]), was simultaneously fit along with our input 
stellar populations in order to correct for contamination 
by unresolved galaxies and foregound stars in our Leo IV 
CMD. 

Stars with magnitudes 18.0 < r < 24.75 (correspond- 
ing roughly to our 50% completeness limit) and —0.50 < 
g — r < 1.15 were fit. After some experimentation, we 
used a Hess diagram bin size of 0.15 mag in magnitude 
and 0.15 mag in color. We assume a binary fraction of 
0.5 and a Salpeter initial mass function. Because we cor- 
rect our stell ar catalog for Galact ic extinction with the 
dust maps of ISchlegel et al.l (|1998l ). we do not allow the 
mean extinction of our model CMDs to vary. As in the 
rest of this work, we chose to fix the distance modulus 
in the code to (to — M)=20.94 mag, although our results 
are robust with respect to this assumption, as we discuss 
below. 

We show the best-fit StarFISH solution in Figure El 
along with a comparison of the best model CMD with 
that of Leo IV in Figure [101 Leo IV is consistent with 
a single stellar population with an age of ~14 Gyr and 
a [Fe/H]=-2.3, although the error bars and upper lim- 
its indicate that there is latitude for both a small, young 
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stellar population and a mix of metallicities at older ages. 
Thus, despite the visual impression of a 'thick' giant 
branch, our analysis does not require a metallicity spread 
to match the oberved CMD. Thorough spectroscopy of 
all stars in the red giant region will be required to deter- 
mine Leo IV membership and to properly quantify the 
metallicity spread. This result is robust to small changes 
in the distance modulus. If we alter the input distance 
modulus from (m - Af)=20.94 to (m - A/)=20.84, than 
a larger fraction of the best-fitting SFH comes from the 
[Fe/H]=-1.3 bin, although the [Fe/H]=-2.3 bin still dom- 
inates. Likewise, a distance modulus of (m — M)=21.04 
yields an old, [Fe/H]=-2.3 stellar population with even 
less from the [Fe/H]=-1.7 and -1.3 bins. 

The model and observed CMDs match remarkably well 
(Figure I10[) given that there are known mismatches be- 
tween the theoretical is ochrones and em pirical, single 
population CMDs (e.g., iGirardi et al.l [20041 ). Also, the 
available models do not span the m etallicity ran g e that 
is apparent in the new MW dwarfs; iKirbv et al.l |2008) 
found ([Fe/H]) = —2.58 with an intrinsic scatter of 
a [Fe/H] = 0.75 in Leo IV, while the Girardi isochrone 
set reaches down to [Fe/H] =-2.3. There is a slight mis- 
match in the BHB between model and data, with the 
best-fitting model CMD having a BHB which is ^0.1- 
0.2 mag brighter than that observed (this is a factor of 
~2 larger than the small BHB magnitude correction we 
made in § 13. 3j) . This could be due to a slight error in our 
assumed distance modulus or a true stellar population 
that is even more metal poor than the most metal-poor 
model in the Girardi isochrone s et, which we know to be 
the case from lKirbv et al.l ((2008). Nonetheless, our basic 
finding that Leo IV is predominantly old (^14 Gyr) and 
metal poor is unaffected. 

4.2. Evidence for Young Stars 

It is well known that dwarf spheroidals harbor a popu- 
lation of blue stragglers - a hot, blue extension of ob- 
jects which lie al ong the normal main sequence (e.g., 
iMateo et al.lll995t ). Because the densities in dwarfs do 
not reach that necessary to produce collisional binaries 
(as they can in the cores of globular clusters) , it is likely 
that they are primordial binary systems, coeval with the 
bulk of stars in the dwarf. Unfortunately, their posi- 
tion along the main sequence makes it difficult to dis- 
entangle blue straggler stars from young main sequence 
stars in the MW's dwarf spheroidals, which has been a 
continuous source of am biguity (e.g.. lMateo et al.l 119951 : 
iMapelli et aHI^O^l [20091 among many others). It is very 
difficult to exclude the possibility that some of the blue 
plume stars in any given dwarf are actually young main 
sequence objects. 

We now articulate two arguments in support of the hy- 
pothesis that at least some of the stars in the blue plume 
of Leo IV are young. First, the blue plume stars appear 
to be segregated within the body of Leo IV. We plot the 
position of the 'high probability' blue plume stars with 
low background/foreground contamination, as identified 
within the dotted box in Figure [TJ onto our smoothed 
map of Leo IV (Figure ITT]) . All of the selected high- 
probability blue plume stars within the body of Leo IV 
are on one side. If Leo IV is assumed to be spherically 
symmetric, the chances of all seven being on the same 
half of the galaxy are ~ (1/2 7 ) or less than 1 percent. 



This segregation would be difficult to explain if all of 
these objects were blue stragglers, which would presum- 
ably have the same distribution as the galaxy as a whole. 
We understand that this a posteriori argument is insuffi- 
cient on its own, but it should be considered in the con- 
text of the high normalized blue plume fraction, which 
we now discuss. 

Our second argument in favor of a young population of 
stars in Leo IV stems from the high blue plume frequency 
no rmalized by the BHB star cou nts, following the work 
of iMomanv et all ([2007D . Briefly. IMomanv et all (|2007l ) 
sought to explore the ambiguity between young main se- 
quence stars and genuine blue straggler stars in a sample 
of MW dwarf galaxies by calculating the number of to- 
tal blue plume objects with respect to a reference stellar 
population - the BHB. The basic result of their work, 
whose data was kindly provided by Y. Momany, was that 
those dwarf spheroidals that do not have a true, young 
stellar component have a blue plume fraction that follows 
a relatively well defined My vs Io^Nbp/Nbhb) anti- 
correlation, where the blue plume consists of only blue 
straggler stars (Figure [T2]) . The physical origin of this 
anti-c orrelation is unclear (see, however, iDavies et al.1 
12004 for a plausible explanation for the anti-correlation 
in globul ar clusters), although it is seen in both open 
clusters jd e Marc hiet al.l [2006) and globular clusters 
(jPiotto et al.1 [2004) . Carina, the one dwarf galaxy in 
their sample which does ha ve a known, recent bout 
of star formation ( ~l-3 Gyr; iHurlev-Keller et all 119981 : 
iMonelli et all2003l ). was shown to have an excess of blue 
plume stars with respect to the aforementioned relation 
(shown as a lower limit in Figure [T2l due to difficulties in 
accounting for blue stragglers associated with the older 
and fainter main sequence turn-off in that system), point- 
ing to the fact that Carina's blue plume was populated 
with young main sequence stars in addition to a standard 
blue st raggler population. More recently, iMartin et al.1 
( 2008a) also found a high blue plume frequency in Canes 
Venatic i I, which we also show in Figure[T2J IMartin et al.l 
(2008a|) used the combination of spatial segregation and 
blue plume frequency to argue that Canes Venatici I har- 
bors a young stellar component. 

We take the ratio of blue plume to BHB stars as iden- 
tified in Figure [TJ by the dashed and solid regions, re- 
spectively, utilizing those stars within the half light ra- 
dius of Leo IV and making background and complet- 
ness corrections. We calculate a blue plume frequency 
of log(A BP /A BHB )=0.56±0.13 for Leo IV and plot this 
ratio along with those of the MW dwarf spheroidals just 
discussed in Figure[12j As can be seen, Leo IV lies off the 
standard My vs log(A t bp / 'A t bhb) anti-correlation just as 
Canes Venatici I and Carina do. Leo IV lies 2 — a a way 
from the linear relation fit by iMomany et al.l (|2007f ) for 
the non-star forming dwarfs. 

Taken together, the segregation of blue plume stars 
along with the high blue plume to BHB fraction points 
to at least some of the stars being young main sequence 
objects. We know that at least one of these blue plume 
stars is a blue straggler, due to the discovery of o ne SX 
Phoenicis variable in Leo IV (|Moretti et al.H2009l) . This 
is only the third of the new dwarf spheroidals (excluding 
Leo T, which appears to be a transition object) that har- 
bor a young stellar population - the others being the pre- 
viously discussed Canes Venatici I ([Martin et al.l [2008a') 
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and Ursa Major II (|de Jong et al.ll2008bl) . We note that 
our evidence for young stars is similar to that presented 
for Canes Venatici I, for which it was also found that the 
blue plume popul ation is segregated a nd blue straggler 
frequency is high ([Martin et al.ll2008al ). We measure the 
luminosity of the young stars in § 14 31 and discuss the 
implications of Leo IV's young stellar component in § [SJ 

4.3. Absolute Magnitude 

As pointed out by iMartin et all (|2008bl ). measuring 
the total magnitudes of the new MW satellites is diffi- 
cult due to the small number of stars at detectable lev- 
els. We account for this 'CMD shot noise' by mimick- 
ing our measurem ent of the total magnitude of Hercules 
(ISand et a l. 2009), which borrowed heavily from the orig- 
inal IMartin et all (|2008bD analysis. 

We take the best SFH solution presented in § 14. 1[ and 
create a well-populated CMD (of ~200000 stars) incorpo- 
rating our photometric completeness and uncertainties, 
using the repop program within the StarFISH software 
suite. We drew one thousand random realizations of the 
Leo IV CMD with an identical number of stars as we 
found for our exponential profile fit, and determined the 
'observed' magnitude of each realization above our 90% 
completeness limit. We then accounted for those stars 
below our 90% completeness leve l by using luminosit y 
function corrections derived from iGirardi et ahl (|2004l ). 
using an isochrone with a 15 Gyr age and [Fe/H]=-2.3. 
We take the median value of our one thousand random 
realizations as our measure of the absolute magnitude 
and its standard deviation as our uncertainty (Table [2]). 
To convert from M r magnitudes to My magnitudes we 
use V - r = 0.16 (TWalsh et al.ll2008l) . 

We find My = — 5.5 ± 0.3 and a central surface bright- 
ness, assuming our exponential profile fit, of /Xo,V = 
27.2±0.3. Both our total absolute magnitude and central 
surface brightness measur ement agrees t o withi n 1 — a 
with the measurement of IMartin et al.l (|2008b| ) , which 
used only SDSS data. 

We use a similar technique for determining the approx- 
imate luminosity of Leo IV's young stellar population 
only. We draw an equivalent number stars from our well- 
populated StarFISH CMD as in Leo IV within our high 
probability blue plume box (see Figure [Hand § I4.2|I . and 
corrected for the luminosity of stars outside this region 
using a luminosity function derived from an isochrone 
with an age of 1.6 Gyr and [Fe/H]=-1.3 ([Girardi et al.l 
2004). We again draw 1000 random realizations to de- 
termine our uncertainties. We find that the young stellar 
population has My — —2.1 ± 0.5, or roughly ^5% of the 
satellite's total luminosity or ~2% of its stellar mass. 
This magnitude and the resulting fraction of the young 
stellar populations luminosity should be taken with cau- 
tion given our assumptions and the small number of stars 
involved. 

5. DISCUSSION & CONCLUSIONS 

In this work we have presented deep imaging of the 
Leo IV MW satellite with Megacam on the MMT and 
study this galaxy's structure and SFH. In particular, we 
assess reports in the literature concerning both its stellar 
population and its possible association with the nearby 
satellite, Leo V. 

Leo IV's SFH is dominated by an old (> 12 Gyr), 



metal poor ([Fe/H] < —2.0) stellar population, although 
we uncover evidence for a young sprinkling of star forma- 
tion 1-2 Gyrs ago. Our best-fit StarFISH results indicate 
that a single metal poor population dominates, although 
the data is also compatible with a spread in metallici- 
ties. The old population is consistent with the emerging 
picture that the faintest MW satellites are 'reionization 
fossil s' (e.g. iRicotti fc GnedirJl2005t iGnedin fc Kravtsovl 
2006), who formed their stars before reionization and 
then lost most of their baryons due to photoevaporation. 
The apparent sprinkling of young stars begs the question 
of what has enabled Leo IV to continue forming stars at a 
low level. There is no sign of HI in Leo IV, with an upper 
limit of 609 M Q (|Grcevich fc Putmanll2009l) . although we 
note that this limit is still a factor of ~2 larger than the 
stellar mass associated with the young stellar population 
studied in § ET2"1 

One possible mechanism for late gas accretion, and 
subsequent star formation, amo ng the f a int M W satel- 
lites was recently discussed by iRicottil <|2009f ) to help 
explain the complex SFH and gas content o f Leo T 
(|de Jong et al.ll2008al: IRvan- Weber et al.ll2008[ ). In this 
scenario, the smallest halos stop accreting gas after reion- 
ization as expected, but as their temperature decreases 
and dark matter concentration increases with decreasing 
rcdshift they are again able to accrete gas from the inter- 
galactic medium at late times, assuming they themselves 
are not accreted by their parent halo until z < 1 — 2. 
This can lead to a bimodality in the SFH of the satellite, 
with both a > 12 Gyr population and one that is < 10 
Gyr, as we see in Leo IV. One stringent requirement of 
this model is that the satellite can not have been ac- 
creted by its host halo until z < 1 — 2 (and thus not 
exposed to tidal stirring and ram pressure stripping un- 
til late times, allowing the satellite to retain its newly 
accreted gas). Future proper motion studies will be able 
to test if Leo IV is compatible with this late gas accretion 
model. Another prediction of this model is the possible 
existence of gas-rich minihaloes that never formed stars, 
but could serve as fuel for star formation if they encoun- 
tered one of the luminous dwarfs. More detailed study 
of this late gas accretion mechanism will be necessary 
to understand the possible diversity of SFHs in the faint 
MW satellites. 

Additionally, we note that if the apparent segregation 
of young stars in Leo IV is real, then it is not an isolated 
case am ong the MW satellite s. As has been mentioned 
in § [4721 IMartin et all (|2008af ) noted that Canes Venatici 
I has a compact star forming region clearly offset from 
the galaxy as a whole, with an age of 1-2 Gyrs, sim- 
ilar to Leo IV. Additionally, Fornax has several com- 
pact clumps and shells that house young stellar popula- 
tions roughly ~1.4 Gy rs old ([Coleman et al.ll2004l 120051: 
lOlszewski et al.ll2006| ). It has been suggested that these 
could have been the result of a collision between For- 
nax and a low-mass ha lo, which was possibly gas-rich. 
iPenarrubia et all (|2009h investigated the disruption of 
star clusters in triaxial, dwarf-sized halos and found that 
segregated structures can persist depending on the or- 
bital properties of the cluster, providing yet another vi- 
able mechanism. More work is needed to distinguish be- 
tween all of the above scenarios and to properly model 
the emerging diversity of SFHs among the new, faint MW 
satellites. 
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Structurally, Leo IV appears to be very round, with 
e < 0.23 (at the 68% confidence limit) and a half light 
radius (~ 130 pc) which is typical of the new MW satel- 
lites. An exhaustive search for signs of extended struc- 
ture in the plane of the sky has ruled out any associated 
streams with surface brightnesses of \i r < 29.6. The ex- 
tent of Leo IV along the line of sight is less than ~15 kpc, 
a limit that will be difficult to improve upon given the 
inherent limitations of using the spread in BHB magni- 
tudes to measure depth. We find no evidence for struc- 
tural anomalies or tidal disruption in Leo IV. We do not 
have the combination of image depth and area necessary 
to confirm the stellar bridge, with a surface brightness 
of 32 mag ar csec" 2 , recently repo rted in between Leo IV 
and Leo V ()de Jong et alj 12010). Indeed, Leo V is al- 
most certainly disrupting, as discussed by iWalker et al.1 
(2009), due to the presence of two member stars ~13 
arcminutes (> 1377,) from Leo V's center along the line 
connecting the putative Leo IV/Leo V system. The na- 
ture of Leo V is still very ambiguous, with the kinematic 
data being consistent with it being dark matter free - 
suggesti ng that perhaps L eo V is an evaporating star 
cluster ([Walker et al.ll2009| ). It is thus critical to obtain 
yet deeper data on these two systems and the region sep- 
arating them to uncover their true nature. 

The probable detection of a small population of young 
stars illustrates once again that it is crucial to obtain 
deep and wide field follow up for all of the newly detected 
MW satellites. Every new object has a surprise or two 
in store upon closer inspection. 
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7. APPENDIX 

In this Appendix, we determine how well the maxi- 
mum likelihood technique presented in ij |3.1l can measure 
the structural parameters of a dwarf with a comparable 
number of stars as Leo IV. We create mock models of 
Leo IV-like systems having an exponential profile with 
77,,=3.0 arcmin and 400 stars in the input catalog (all on 
a footprint with the same area as our MMT pointing), 
while systematically varying the ellipticity and position 
angle. A uniform background of 4.5 stars per square ar- 
cminute is randomly scattered across the field to mimic 
the actual observations. We utilize the same algorithm 
as in § 13.11 with 1000 bootstrap resamples to determine 
our uncertainties. 

We can recover the ellipticity of our mock dwarf galax- 



ies remarkably well, as can be seen in Figure H3l Here we 
show our results on the recovered ellipticity as a function 
of input model ellipticity, between 6i npu t =0.05 and 0.85. 
The data points are the median ellipticity found for the 
one thousand bootstrap resamples for each model, and 
the error bars encompass 68% of the resamples around 
that median. Note that for models with an input el- 
lipticity of e < 0.25 it is difficult for our algorithm to 
converge on the correct ellipticity value. In this regime, 
we present the measured ellipticity as an upper limit (see 
§ 13.11) . In this case, we systematically overpredict the el- 
lipticity with large error bars, and thus only quote upper 
limits. Larger values of the ellipticity are well measured. 

We also do a good job of measuring the position an- 
gle (as long as the ellipticity is high enough) and the 
half light radius of our mock Leo IV-like systems, which 
we illustrate in a series of examples in Figure 1141 From 
the figure, one can see the gradual improvement in the 
measurement of the position angle as one goes from an el- 
lipticity of 0.1 to 0.35, while the half light radius remains 
relatively well measured throughout. The slight system- 
atic offset seen in the bottom left panel of Figure [14] can 
be explained due to the difficulty in recovering the true 
ellipticity for ei npu t < 0.25 systems. If one slightly over- 
estimates the true input ellipticity of the data, this leads 
to a slight overestimation of the half-light radius, a de- 
generac y that can be seen in Figure 9 of iMartin et al.1 
(2008b]). The take away message is that one must treat 
any 'measurement' of the position angle with extreme 
caution for ellipticity values of e < 0.25 as we do in § 13.11 

In the future, we intend to present a more extensive 
series of tests of our maximum likelihood code in order 
to understand how star number and imaging field of view 
have on the estimation of structural parameters for the 
new MW satellites. 
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Fig. 1. — Left - The CMD of Leo IV within 7"h=2.85 arcmin. We note several regions marking off particular stellar populations. The 
dash-lined region encloses the total blue plume population, although any given star may be a foreground contaminant or unresolved galaxy. 
The inset dotted box is our high probability blue plume box, with relatively little galactic contamination. Finally, we consider those stellar 
objects w ithi n the solid-lined region to be BHB stars. We plot the spatial distribution of the high probability blue plume stars and the BHBs 
in Figure im Right - A bac kground subtracted Hess diagram of the same, central half light radius region of Leo IV. Overplotted are several 
theoretical isochrones from Girardi et al. (2004). The solid and dashed ridgelines are of a 14 Gyr stellar population with [Fc/H]=-2.3, -1.7, 
respectively. The dotted line is a ~1.6 Gyr isochrone with [Fe/H]=-1.3, and roughly matches the overdensity of blue plume stars that are 
evident. 
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Fig. 2. — An illustration of the CMD selection region used for performing our parameterized structural analysis of Leo IV. The figure 
is a background-subtracted Hess diagram of the central 2.5 arcminutes of Leo IV and the central, red fiducial ridge line is that of M92 
transformed to a distance modulus of (m — M) = 20.94. The two bordering ridge lines bound the selection region and arc a minimum of 
0.1 magnitudes in the color direction away from the M92 ridge line. The selection region expands to the size of the typical g — r color 
uncertainty as a function of r mag when this exceeds 0.1 mag, as determined via our artificial star tests. The selected region is an excellent 
match to the apparent ridgeline of Leo IV, as can be seen in the figure. 
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Fig. 3. — Stellar profile of Leo IV. The data points are the binned star counts for all stars within the CMD selection box shown in Figure[2] 
The plotted lines show the best fit one-dimensional exponential, Plummer and King profiles. The dotted line shows the background surface 
density determined for our exponential profile fit. Note that in deriving these best fits, we are not fitting to the binned data, but directly 
to the stellar distribution. 
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Fig. 4. — Smoothed contour plots of Leo IV, showing the 3, 4, 5, 6, 7, 10 and 15 a levels. The left panel has been smoothed by 0.5 
aremin, the middle panel by 1.0 aremin, and the right panel by 1.5 aremin. The arrow in the middle panel indicates the direction to Leo V. 
Note that there is no apparent structure along this direction. The 3 a contour in our 0.5 arcminute smoothed map corresponds roughly to 
a surface brightness of ~30.0 mag per arcsccond 2 . 
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Fig. 5. — Smoothed contour plots of nine random realizations of Leo IV stars, where we have reassigned star positions randomly across 
the Megacam field of view. The contours show the 3, 4 and 5 a levels. The plots show that compact 3 and 4 a overdensities are relatively 
common. 
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Fig. 6. — Smoothed contour plots of 9 random bootstrap resamples of Leo IV stars (with the 0.5 arcminute Gaussian), showing how 
morphology can vary with each iteration. The contours show the 3, 4, 5, 6, 7, 10 and 15 a levels. While in some iterations the 'tendrils' 
seen in Figure [3] are still visible, they are not a ubiquitous feature. Future, deeper imaging is necessary to confirm their reality. 
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FlG. 7. — An illustration of our technique for implanting fake Leo IV 'nuggets' into our field, using an artificial stellar population drawn 
from that determined in our star formation history analysis I4.1[l and our artificial star tests (§E). Left: We show our smoothed map 
of Leo IV after implanting a 35 star 'nugget' at (+5', +5') with respect to the center of Leo IV, distributed as an exponential profile with 
a half-light radius of 1 arcminute (see Table [2J. This nugget results in a ~6.4-<r overdensity at that position. Right: A Hess diagram of 
the 35 star nugget, along with a M92 isochrone shifted to (m — M) = 20.94. This Hess diagram shows the raw star counts in the circular 
aperture shown in the left panel of this figure, along with an equal area background region, and the difference between the two. Note that 
the residual CMD shows several stars that would be identified with Leo IV's red giant branch, which would allow us to say with confidence 
that such an overdensity is likely associated with Leo IV. 
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Fig. 8. — An illustration of our technique for implanting fake Leo IV 'streams' into our field. Top: We show our smoothed maps (using 
a 1 arcminute Gaussian) of Leo IV after implanting streams with different numbers of stars (see Table |3). Note that the 200 and 300 star 
streams are clearly detected. Bottom: The resulting Hess diagrams for our 200 and 300 star scenarios. Note that each has many obvious 
red giant and BHB members, while the 300 star stream has the beginnings of the main sequence. We show the region we used for extracting 
the 'stream' photometry in our maps in the top row of this figure. 
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Fig. 9. — SFH solution from the StarFISH fit. With the data in hand, Leo IV is consistent with a single, old (~14 Gyr) stellar population 
with [Fe/H]~-2.3 - although there is some weak evidence for a spread in metallicity in this old stellar component. We also can not rule 
out a small, young population of stars with our StarFISH analysis. Error bars with no accompanying histogram are upper limits. 
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Fig. 10. — Comparison of the data with our best StarFISH model fit. Pixel bins are 0.15 mag along the color and magnitude axis (except 
for panel (d), where the bins are 0.05 mag), as used in our StarFISH fits. The gray scale for each Hess diagram is in units of stars per pixel, 
except for panel (f), which has been scaled by the uncertainty associated with each pixel, (a) The observed Hess diagram of Leo IV within 
the half light radius, (b) The background Hess diagram, which was fit along with the theoretical isochrones. (c) The background-subtracted 
Hess diagram, (d) The best model CMD derived by StarFISH (with no background component), in bins of 0.05 mag in order to show the 
details of the model, (e) Residuals after subtracting a random realization of the StarFISH model from the data, (f) The residual from (e), 
scaled by the expected Poisson scatter in the bin. The most significant residuals are associated with the mismatch of the BHB between 
model and data. 
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Fig. 11. — The position of our high-probability blue plume (squares) and candidate BHB (diamonds) stars overplotted onto the smoothed 
map of Leo IV. Note that the probable blue plume stars are segregated within the body of Leo IV, as would be expected if they represented 
a young stellar population rather than blue straggler stars. As expected, there are few of these objects outside the main body of Leo IV. 
The BHB star population is more uniformly distributed across Leo IV. 
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Fig. 12. — The frequency of blue plume stars, normalized by BHB stars, for the M W dwarf spheroidals. The solid, circular points are the 
blue plume frequency points for the MW dwarfs derived by Momanv ct al. (2007]). Note the anti-correlation between satellite brightness 
versus blue p lume frequency. An outlier from that work is Carina, which likely harbors a young stellar population. The Canes Venatici I 
point is from Martin ct al. ( 2008a:) who reported a young stellar population in that system due to both the high blue plum e frequency and 
the se gregation of this population within the dwarf. Li kew ise, the blue plume frequency of Leo IV lies off the relation of Momanv ct al. 
(2007). That, along with the segregation seen in Figure fTTI indicates that Leo IV harbors a small, young stellar population. 
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Input ellipticity 

Fig. 13. — Comparison of the input ellipcity in our mock Leo IV dwarf galaxy models versus the output ellipticity found by our maximum 
likelihood algorithm. The points represent the median value of the 1000 bootstrap resamplcs, and the error bars represent the central 
68.3% of the distribution around the median. The dashed line is the one-to-one correspondence between input and output. For input 
ellipticity values of ei npu t < 0.25 we systematically overpredict the ellipticity with large error bars and thus only quote upper limits. 



Leo IV SFH and Structure 



23 







Position Angla (deg) 



n 









Fig. 14. — An illustration of how the maximum likelihood code for measuring parameterized structure performs for an input position 
angle of 35 degrees while varying the input ellipticity between 0.1 and 0.35. The output histograms are the results of the one thousand 
bootstrap rcsamples for the labeled scenarios, while the solid vertical line is the input value of the model, in this case PA = 35.0 degrees 
(top row) and = 3.0 arcmin (bottom row). Top row - Recovered position angle for various mock Leo IV models with e between 0.1 
and 0.35. Note how the example with an input ellipticity of e=0.10 does not have a clearly measured position angle, and the peak of the 
distribution docs not agree with the model's position angle. It is not until an ellipticity of e ~0.35 where the model's position angle is 
well measured and accurate. Bottom - Recovered half light radius while varying the input ellipticity. Although the low ellipticity case of 
e = 0.1 seems to be slightly offset from the input = 3.0 arcmin model, the recovered half light radius is still within l-cr of the input 
value. In all other cases, the half light radius is measured with high precision. 
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TABLE 1 
Leo IV Photometry 



Star No. 


a 


<5 


9 






r 


<5r 


A r 


SDSS or MMT 




(deg J2000.0) 


(deg J2000.0) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 







173.23115 


-0.54635 


18.23 


0.02 


0.09 


17.93 


0.01 


0.07 


SDSS 


1 


173.24960 


-0.50828 


17.16 


0.01 


0.09 


16.61 


0.01 


0.06 


SDSS 


2 


173.20208 


-0.54312 


17.03 


0.02 


0.09 


16.61 


0.01 


0.07 


SDSS 


3 


173.22927 


-0.49215 


18.66 


0.02 


0.09 


17.20 


0.01 


0.06 


SDSS 


4 


173.25937 


-0.50253 


17.87 


0.01 


0.09 


17.47 


0.01 


0.06 


SDSS 


5 


173.19711 


-0.49921 


16.58 


0.02 


0.09 


15.84 


0.01 


0.06 


SDSS 


6 


173.21481 


-0.58547 


16.52 


0.02 


0.09 


15.50 


0.01 


0.07 


SDSS 


7 


173.23169 


-0.47729 


17.00 


0.02 


0.09 


15.81 


0.01 


0.06 


SDSS 


8 


173.23995 


-0.60117 


18.08 


0.02 


0.09 


17.47 


0.01 


0.07 


SDSS 


9 


173.19902 


-0.56584 


19.31 


0.02 


0.09 


17.78 


0.01 


0.07 


SDSS 



See electronic edition for complete data table. 
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TABLE 2 

LEO IV STRUCTURE - PARAMETERIZED FITS 



Parameter 


Measured 


Uncertainty 


Bootstrap median 


My 


-5.5 


0.3 


-5.5 




27.2 


0.3 


27.2 


Exponential Profile 


RA (h m s) 


11:32:56.38 


±14" 


11:32:56.33 


DEC (d m s) 


-00:32:27.25 


±14" 


-00:32:29.24 


(arcmin) 


2.85 


0.64 


3.17 


(PC) 


127.8 


28.8 


142.2 


e 


0.05 


< 0.23 


0.14 


Plummer Profile 


RA (h m s) 


11:32:56.20 


±12" 


11:32:56.18 


DEC (d m s) 


-00:32:27.17 


±10" 


-00:32:28.32 


(arcmin) 


2.86 


0.40 


2.97 


(PC) 


128.1 


18.0 


133.3 


€ 


0.02 


< 0.20 


0.11 


King Profile 


RA (h m s) 


11:32:56.04 


±13" 


11:32:56.08 


DEC (d m s) 


-00:32:29.60 


±10" 


-00:32:29.95 


r c (arcmin) 


1.61 


0.22 


1.64 


(PC) 


72.2 


10.3 


73.4 


rt (arcmin) 


18.55 


4.39 


20.25 


(PC) 


831.0 


196.7 


907.3 


e 


0.03 


< 0.18 


0.09 



a All transverse distances are reported using a (m — M) — 20.94 
distance modulus. 

b Absolute magnitude and central surface brightness are calcu- 
lated using the exponential profile fit. 

c The value in the uncertainty column for e corresponds to the 
68% upper confidence limit, given that our derived e is consistent 
with 0. 
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TABLE 3 

Simulated Leo IV External Structures and 
Detections 



No. of stars M r /io.r M g /no,g Peak a 



Input 'Nuggets' 



25 


-1.0 


29.7 


-0.6 


30.1 


4.0 


35 


-2.8 


27.9 


-2.5 


28.3 


6.4 


50 


-3.2 


27.5 


-3.0 


27.7 


9.7 


Input 'Streams' 












100 


-3.5 


30.5 


-3.2 


30.8 


1.7 


200 


-4.4 


29.6 


-4.3 


29.8 


4.3 


300 


-5.4 


28.7 


-5.0 


29.0 


5.1 





a All nuggets have an exponential profile with half 
light radius of 1 arcminutes. 

b All streams have a constant density in the dec- 
lination direction, with a gaussian density profile 
in the right ascension direction, with a— 1.5 ar- 
cmin. 

c We quote the peak a for our artificial streams at 
a position of (+0.0, +8.0) arcminutes with respect 
to the center of Leo IV 



